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ABSTRACT 

We analyzed 123 thermonuclear (type-I) X-ray bursts observed by the Rossi X-ray Timing Explorer 
from the low-mass X-ray binary 4U 1636—536. All but two of the 40 radius-exansion bursts in this 
■ sample reached peak fluxes which were normally distributed about a mean of 6.4 x 10~^ ergs cm~^ s~^ , 

' with a standard deviation of 7.6%. The remaining two radius-expansion bursts reached peak fluxes a 

factor of 1.69 ± 0.13 lower than this mean value; as a consequence, the overall variation in the peak 
flux of the radius-expansion bursts was a factor of « 2. 
^ ' This variation is comparable to the range of the Eddington limit between material with solar H- 

O , fraction {X — 0.7) and pure He. Such a variation may arise if, for the bright radius-expansion bursts, 

most of the accreted H is eliminated either by steady hot CNO burning or expelled in a radiatively- 
drivcn wind. However, steady burning cannot exhaust the accreted H for solar composition material 
within the typical « 2 hr burst recurrence time, nor can it result in sufficient elemental stratification to 
allow selective ejection of the H only. An additional stratification mechanism appears to be required to 
separate the accreted elements and thus allow preferential ejection of the hydrogen. We also observed 
non-radius expansion bursts that exceeded the peak flux of the faintest radius expansion bursts. For 
' these bursts the accreted hydrogen must have been partly ejected or eliminated, but the burst flux 

I did not subsequently reach the (higher) Eddington limit for the underlying He-rich material. 

, We found no evidence for a gap in the peak flux distribution between the radius-expansion and 

non-radius expansion bursts, previously observed in smaller samples. Assuming that the faint radius- 
expansion bursts reached the Eddington limit for H-rich material {X « 0.7), and the brighter bursts 
lO i the limit for pure He {X — 0), we estimate the distance to 4U 1636—536 (for a canonical neutron star 

O ■ with Mns = 1.4Mq, i?NS = 10 km) to be 6.0 ± 0.5 kpc, or for Mns = ^Mq at most 7.1 kpc. 

, Subject headings: stars: neutron — X-rays: bursts — stars: individual (4U 1636—536) — stars: 

pi i' distances 

o 

|H I L INTRODUCTION ignition, X (as well as the temperature of the fuel layer) 

^ ■ Thermonuclear (type I) X-ray bursts are caused by un- modified by steady hot-CNO H-burning between the 

?^ ■ stable burning of accreted matter on the surface of neu- H'n'''*'' (e-g-| Fni™"to et al. || 198l D . Such bursts have been 

tron stars in l ow-mass X-ray binary (LMXB) systems observed to date from more than 70 sources (for a recent 

^ (seeEiiEiEal. 1993; Bildsten 1998, for reviews). Typ- catalog see| m t Zand et al | 2004D . , , ^ , 

<N ■ ical burst profiles exibit short rise times between < 1 to energy from the burst is released suSiciently 

^ ^ 10 s, and decay time scales from 10 to ~ 100 s. Model fits ^^pidly the flux niay exceed the Eddington limit which 

■ ■ ' using a blackbody continuum to X-ray spectra during the ^""^ ^ suSiciently distant observer is ( |Lewm et al. || 1993D 

bursts provide evidence for an initial rise in color temper- 87rGrnpMNsc[l -t- (aTll;)"-®^] / 2GM^s\^^'^ 

ature Tbb, followed by a more gradual decrease back to ^Edd^oo = - ^ ^ -^^ I 1 I 

persistent levels. This is naturally interpreted as heating 

resulting from thermonuclear ignition of surface fuel, fol- — 3 5 x 10'^^ I ] 1 + (qtJc. 
lowed by cooling of the ashes once the available fuel is ex- ' \1.4Afo/ 1 + X 
hausted. The inferred blackbody radius is around 10 km, \ 1/2 
consistent with expectations for a wide range of neutron- ^ f 1 ^^-^"^ns \ ergss"^ 
star equations of state. The time to achieve burst igni- V Rc^ J 
tion depends primarily on the accretion rate and the H- where Mns is the mass of the neutron star, is the ef- 
fraction, Xo, in the accreted material; the composition at fective temperature of the atmosphere, ax is a coefficient 

parametrizing the T-dependence of the electron scatter- 
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sufRcient to lift the outer layers of the atmosphere above 
the neutron star surface. Thermonuclear bursts exhibit- 
ing photospheric radius-expansion (PRE) are thus im- 
portant because the peak flux can be estimated based on 
the gravitational redshift and atmospheric composition. 
If the emission is isotropic, such bursts represent a "stan- 
dard candle", which can in princi ple allow estimates o f 
the distance to the source (e.g. see lKuulkers et al1l2003H 
or th e compactness of the neutron star fe.g. TDamener^Ll 
[1993). 

In this series of papers, we investigate empirically the 
assumption that the peak flux of PRE bursts is constant 
for each burst source. In the first paper, we found ev- 
idence for significant variation in the peak flux of PRE 
bursts observed with the Rossi X-ray T iming Explorer 
(RXTE) from 4U 1728-34 (Galloway et al . 200_3a, here- 
after Paper A). The peak burst fluxes appeared to vary 
steadily on a timescale of a few tens of days, which was 
similar to the timescale at which the persistent X-ray 
flux was modulated. The « 10% rms variation in peak 
burst flux was attributed to varying degrees of reflection 
from a precessing accretion disk, with the burst emission 
inferred to be intrinsically isotropic. Here we present a 
study of the variation of the peak fluxes of PRE bursts 
from 4U 1636—536, through analysis of the largest sam- 
ple to date, gathered from the available public data from 
observations by RXTE. 

4U 1636-536 {I = 332?9, h = -4?8) is a well-studied 
LMXB, consisting of a neutron star i n a 3.8 h orbit with 
an 18 th magnitude star, V801 Ara l)van Paradiis et all 
119901 see also Giles et al. 20^o3). The X-ray source 
exhibits a variety of rapid tim e variability, includin g 
kHz quasi-periodic oscillations ll Wiinands et alJ 119971 ^ , 
X-ray bursts, and burst oscillations at 579.3 Hz 
ijStrohmaver et aDll998alHl . Previous analyses of small 
numbers of X-ray bursts observed from this source by 
various satellites revealed that their peak fluxes appeared 
to be distributed bimodally in the most part, with the 
PRE bursts reaching a peak flux a f actor of 1.7 higher 
than the brightest non-PRE burst fin oiie et al.l 119841: 
^ewin ct al. 1987). The measured peak fluxes of the PRE 
bursts were generally found to be consistent from burst to 
burst fe.g.. lEbisuzakilll98';1) . The properties of the PRE 
bursts, along with the claimed detectio n of redshifted ab- 
sorption features in the burst spectra ijWaki et alJll984|) 
have led to distance est imates between 6-7 kpc (e.g. 
IChristian 8^317991. 

2. RXTE OBSERVATIONS OF 4U 1636-536 

The High-Energy Astrophysics Science 
Archive Research Centre (HEASARC; 

\protect ittp://heasarc. gsfc.nasa.gov) contains public 
RXTE observations of 4U 1636—536 dating from shortly 
after the launch of the satellite on 1995 December 
30. We extracted all the available public data, which 
at the time of writing includes observations up to 
2004 March 26. We principally used data from the 
Proportional Counter Array (PCA; I.Tahoda et alJll99^ 
aboard RXTE, which consists of five identical gas-filled 
proportional counter units (PCUs) with a total effective 
area of ~ 6000 cm^ and sensitivity to X-ray photons 
in the 2-60 keV range. We have developed a pipeline 
processing system to identify and download newly public 
data from around 70 known bursting sources, including 



4U 1636—536. We generated 1-s binned lightcurves from 
each observation, and then identified highly significant 
single-bin deviations from the mean count rate as burst 
candidates. Each such candidate was visually inspected 
to distinguish from other possible sources of abrupt 
count rate variation, including PCUs being turned on or 
off, or PCU breakdowns. 

Once located, high time- and spectral resolution PCA 
data (where available) covering each burst were pro- 
cessed to obtain full-range spectra within intervals of 
0.25-4 s (with the integration time increasing as the burst 
count rate decays). To take into account gradual varia- 
tions in the PCA gain we generated a response matrix for 
each burst using pcarsp version 10. l'', which is part of 
LHEASOFT release 5.3 (2003 November 17). A persistent 
emission spectrum extracted from a (typically) 16 s in- 
terval prior to the burst was used as the background; this 
approach is well-established as a standard pro cedure in 
X-ray b urst analysis (e.g. , Kuulk ers et al.l200l . although 
see also Ivan Paradiis fc Lewi n 1986*). We estimated the 
persistent flux, i^per, at the time of each burst by fitting 
the background-subtracted spectrum averaged over the 
entire observation (excluding the bursts) with a model 
consisting of an absorbed blackbody and power law. 

We fitted each time-resolved burst spectrum with a 
blackbody model multiplied by a low-energy cutoff rep- 
resenting interstellar absorption with fixed abundances. 
The initial fitting was performed with the absorption col- 
umn density nn free to vary; subsequently, it was fixed 
at the mean value measured over the entire burst for the 
final results. The fiux at the peak of the burst did not 
vary significantly as a function of njj. The bolometric 
fiux at each timestep ti was calculated according to 



= 1.0763 X 10"" r^b,»-K'bb.* ergscm-^s"^ (2) 

where Tbb is the color temperature, iiTbb = 
(^bb,km/c^iokpc)^ is the blackbody normalisation, 
with i?bb,km the apparent radius of the neutron star 
in km for a distance of diokpc = d/(10 kpc). As a 
working definition, we considered that radius expansion 
occurred when 1) the blackbody normalization /Cbb 
reached a (local) maximum close to the time of peak 
fiux; 2) lower values of ii'bb were measured following 
the maximum, with the decrease significant to a level of 
4(7 or more; and 3) there was evidence of a significant 
(local) decrease in the fitted temperature Tbb at the 
same time as the increase in i^bb (these criteria are 
identical to those used in Paper A). We measured the 
fiuence i?b by integrating numerically over the measured 
values of i^boi.i, extrapolating the derived exponential 
decay curve for the cases where the burst emission lasted 

^ We note that the geometric area of the PCUs was changed 
for this release for improved consistency between PCUs and 
(e.g.) canonical models of calibration sources, particularly the 
Crab pulsar and nebula. These changes have the effect of 
reducing the measured flux compared to analyses using previ- 
ous versions of the response generating tools by 12—14%. S ee 
\protect ittp:/ /Iheawww.gsfc. nasa. gov / ~keith / p ca-calibration_draft .ps| 
(Jahoda et al. 2UU5, ApJ, submitted] tor more details. 
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Fig. 1. — Maximum blackbody radius (for a source distance of 
6.0 kpc) in the rising phase of the burst as a function of the peak 
flux for 78 non-radius expansion and 40 radius expansion bursts 
from 4U 1636—536. Error bars indicate the la uncertainties on 
each measurement; radius-expansion bursts are additio nally indi- 
cated by diamond symbols. Burst (c) is labeled (see i|3.1i . The 
dashed lines delineate the gap in the ^pcak distribution inferred 
by Sujjimoto et al] 119841) . The histograms show the separate peak 
flux distributions of the non-radius expansion (dotted line) and 
radius expansion (solid line) bursts. The largest bin for the two 
distributions contains 18 bursts. 

longer than the high-resolution spectral data (typically 
200 s). 

3. RESULTS 

The currently available public RXTE data include 123 
X-ray bursts from 4U 1636—536. The distribution of 
peak fluxes as a function of the presence or absence of 
radius expansion is shown in Fig. ^ The majority of 
the PRE bursts reached peak fluxes Fp^ak that were dis- 
tributed normally between 5.5-7.4 x 10~^ ergscm~^s~^, 
with a mean of 6.4 x 10^^ ergscm~^s~^ and a stan- 
dard deviation of 7.6%. This variation was significantly 
greater than the typical error on the flux {'^ 3%); a 
calculation indicates that the hypothesis of a constant 
Fpoak for these bursts can be excluded down to a confi- 
dence level of < 10" (equivalent to > 8cr). The PRE 
bursts reached maximum i?bb approximately a factor of 
two larger than the maximum achieved during the rise 
for the non-PRE bursts. The weighted mean maximum 
radii were 14±2 and 6.9±1.0 kmd^Jp^, respectively. The 
smallest maximum _Rbb reached by a burst which showed 
unambiguous indications of PRE was 8.8 kmd^j^p^.. 

We also observed 78 non-PRE bursts, with peak 
fluxes distributed normally between 0.21-6.3 x 
10~^ ergscm"^ s~^. We found several bursts (both PRE 
and non-PRE) which reached peak fluxes between w 3.5- 
6 x 10~^ ergsci r i~^s~ ^, i.e . within the "gap" noted by 
iSugimoto et all l)1984j) and lLewin et alJ l)1987|) . We also 
found two PRE bursts that reached a peak flux of around 
3.7 X 10~* ergscm~^s~^, which is lower than the max- 
imum peak flux of the non-PRE bursts. In the next 
section f i)3.1|l we describe these bursts in detail. 



The 7.6% spread in peak fluxes for the brighter 
(■P'peak^ 5 X 10~^ ergscm~^s~^) PRE bursts is similar 
to that measured in 4U 1728-34 (Paper A). If the vari- 
ation in the peak fluxes of the bright PRE bursts from 
4U 1636—536 arises from the same mechanism, we might 
expect to find evidence for a connection with the persis- 
tent emission, as well as for reprocessing of the burst flux 
(i.e. a correlation between the peak flux and fluence). 
While there are some indications for quasi-periodic vari- 
ation at a period of w 75 d in the long-term All-Sky 
Monitor lightcurve, this periodicity does not appear to 
be reflected by the peak burst fluxes. We found no signif- 
icant peaks in a Lomb-Scargle periodogram of the peak 
fluxes as a function of time and no correlation between 
the peak PRE burst flux and the fluence (as was also 
observed in 4U 1728—34). We also found no evidence 
for modulation of the peak burst fluxes at the 3.8 hr or- 
bital period (usin g the most recent published ephemeris; 
IGiles et a,].ll20?)^ . 

3.1. Radius- expansion at low peak flux 

We found two bursts that exhibited PRE but reached 
peak fluxes signiflcantly below the maximum peak flux 
of the non-PRE bursts. While one of the bursts (on 
1999 Sep 25 20:40:49 UT) exhibited a factor of - 2 in- 
crease in radius and the archetypical flat-topped flux pro- 
file of a PRE burst, the other (on 2000 Jan 22 04:43:48 
UT) had much more modest expansion and an overall 
less compelling case for PRE. The bursts reached peak 
fluxes of 3.76 x 10~^ and 3.60 x 10~^ ergscm~^s~^, re- 
spectively. These peak fluxes were consistent to within 
the errors, and were less than the mean for the remain- 
ing PRE bursts by a factor of 1.69 ± 0.13 (where the 
uncertainty arises primarily from the 7.6% variation in 
the peak fluxes of the latter sample). These bursts rep- 
resent a highly significant deviation from the normally- 
distributed peak fluxes of the remaining PRE bursts, and 
lead to an overall variation of a factor of w 2 in the peak 
fluxes of PRE bursts from 4U 1636-536. 

In Fig. [5]we show the variation of spectral parameters 
throughout a bright PRE burst, a non-PRE burst, and 
one of the fainter PRE bursts. The profiles for the bright 
and faint PRE bursts (a, left panels and c, right panels) 
were similar before and after the PRE episode. However, 
for burst (c) the flux was approximately constant during 
the PRE episode at around 4 x 10~® ergscm~^s~^, i.e. 
at around 60% of that reached at the peak of burst (a) . 
We note that the two faint PRE bursts were observed 
while the persistent 2.5-25 keV flux was unusually high, 
above 6x 10~^ ergs cm~^ s~^. We observed an insufficient 
number of bursts within this flux range to determine con- 
clusively whether the peak fiux distribution differed from 
that at lower persistent fiuxes. 

Both i?bb and Tbb measured at the time of maxi- 
mum fiux of the two faint PRE bursts were, individ- 
ually, within the ranges spanned by the brighter PRE 
bursts. Thus, the lower peak fluxes were not solely due 
to unusually low values of either of these parameters. We 
note that, as with previous observations ijSugimoto et alJ 
'1984!), the radius maximum during the majority of the 
PRE bursts was generally achieved prior to the max- 
imum measured fiux. Thus, the luminosity continued 
to increase throughout the PRE episodes, as has been 
found on other sources fe.g. lGallowav et alJl2003b(l . The 
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two faint PRE bursts were notable exceptions; for those 
bursts, the peak flux coincided with the peak radius. 
Given the relatively wide field of view (~ 1°) of RXTE 
it is conceivable that these fainter PRE bursts actually 
originated from a previously unknown and more distant 
field source. The detection of burst oscillatio ns at 580 Hz 
in the burst on 1999 Se ptember 25 (Fig. iGiles et all 
120021 iMuno et al]l200H) effectively rules out this possi- 
bility, at least for that burst. 

A third burst, observed on 2002 January 15 14:08:16 
UT, showed evidence of a local radius maximum, but at 
an even lower peak flux than burst (c) . During this burst 
the flux rose and decayed gradually, possibly with more 
than one local maxima. The maximum radius reached 
was only somewhat below that of burst (c) and, like 
that burst, was accompanied by a local minimum in Tbb. 
However, the peak flux was below 2 x 10~® ergs cm^^ s^^ , 
i.e. less than 25% the peak flux of burst (a); further- 
more, the possible radius expansion episode occurred 
well before the peak flux was reached, at a flux of only 
around 1.3 x 10~^ ergscm~^s~^. This burst was sub- 
stantially different in character from the other two faint 
PRE bursts. It reached a comparable peak flux and was 
observed at a similar persistent flux level as were four 
double-peaked bursts, also observed by RXTE. These 
bursts (on 2001 September 5 08:15:04 UT, 2001 Octo- 
ber 3 00:22:18 UT, 2002 January 8 12:22:44 UT and 
2002 February 28 23:42:53 UT) exhibited double peaks 
in the bolometric flux, as has been ob served previously 
from this source l)Sztaino et al.lir985^ . In contrast to 
the double-peaked bursts observed previously, in which 
the first peak was consistently higher, 3 of the 4 double 
peaked bursts observed by RXTE reached a higher flux 
during the second peak. The relative fluxes of the two 
peaks varied from w 1-3. Two of the four double-peaked 
bursts exhibited a local maxima in the radius (and a cor- 
responding minimum in temperature) coincident or just 
after the first flux peak, similar to what is generally in- 
terpreted as evidence for PRE. Because of this similarity, 
we conclude that the 2002 January 15 14:08:16 UT burst 
was not a genuine PRE burst, but instead arose from the 
same phenomenon which gives rise to the double-peaked 
bursts. 

4. DISCUSSION 

We have studied the peak fluxes of radius expansion 
bursts from 4U 1636—536 and found the largest range 
yet seen in any LMXB. While the fractional variation 
of peak fluxes exceeded the expected range of the Ed- 
dington limit between atmospheres with H at cosmic 
abundances (^Edd.Hi wit h X ~ 0.7) and with pure He 
(^Edd.He, ^ = 0; e.g. IKuulkers et alJl2C)03f) . the distribu- 
tion of peak fluxes was nearly bimodal, with the mean 
peak flux for the majority of the PRE bursts being a 
factor of 1.69 ± 0.13 higher than the peak fluxes for the 
faint PRE bursts. 

In a previous investigat ion of X-ray bursts from 
4U 1636-536 with Tenma, iSugimoto et'all l|1984H ob- 
served a gap in the distribution of peak fluxes of bursts 
and identified the lower and upper boundaries of this gap 
as the Eddington limits for H-rich LEdd.H and pure He 
fuel iEdd.Hc, respectively. According to their interpeta- 
tion, the increase in radius which was observed during 
the PRE bursts from 4U 1636—536 was accompanied by 
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Fig. 2. — Thermonuclear bursts from 4U 1636—536 showing ra- 
dius expansion over a wide range of peak fluxes. From left to right, 
we show a PRE burst with strong radius expansion, observed on 
1998 Aug 20 03:40:09 UT; a non-PRE burst observed on 2000 Jan 
22 01:46:23 UT; and a faint PRE burst observed on 1999 Septem- 
ber 25 20:40:49 UT. From top to bottom, the panels show for each 
burst the time evolution of the bolometric flux -Fbol (iii units of 
10~^ ergs cm"'^ s~^), the fitted blackbody temperature T^b, the 
blackbody radius i?bb (assuming a source distance of 6.0 kpc), and 
finally the path traced by the burst in Tbb-i^bol space. The cir- 
cles in the lower two panels in each column indicate the time of 
maximal radius expansion, where present. The 3 dashed lines in 
each of the bottom panels represent the expected curves for black- 
body emission from a neutron star with apparent radius 8, 10, and 
15 km (moving left to right, respectively). Error bars represent the 
la uncertainties. 

ejection of an outer, H-rich layer, exposing the underlying 
pure-He layer below. Bursts that never reached iEdd,H 
did not exhibit PRE, while the characteristic flux limit 
for bursts that did exceed £Edd,H switched to iEdd.Hc, 
which is a factor of 1.7 higher. 

In this paper, we reported the observation of two PRE 
bursts that were apparently limited by iEdd,H instead 
of LEdd.Hs. Th i s app ears to confirm the hypothesis of 
ISugimoto et alJ l)1984ft . For those bursts, the effective 
Eddington limit could have remained at LEdd,H if the 
burst energy was insufficient to drive off the outer H- 
rich layer. However, we also found a significant over- 
lap between the peak flux distributions for the PRE and 
non-PRE bursts, without the gap found in earlier sam- 
ples. Bursts with peak fluxes within the range ~ 3.5- 
6 X 10~* ergs cm~^ s""'^ were, however, relatively infre- 
quent. This strongly suggests that the previously ob- 
served gap was merely a consequence of small burst 
samples, specifically 12 fo r the study of ISugimoto et alJ 
(j^84), and 27 for that of iLewin et all l|1987D . with onlv 
3 radius expansion bursts in each sample. 

Four bursts observed by RXTE that peaked between 
the putative iEdd,H and LEdd.Hc values exhibited no evi- 
dence of PRE (an example is shown in Fig. \^). In order 
for these bursts to have exceeded LEdd.H, the majority of 
the accreted H must have been ejected, but subsequently 
the flux must have remained below the effective Edding- 
ton limit for the residual material. The fact that no ob- 
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servational evidence of this ejection is seen suggests that 
the ejected material becomes transparent on a timescale 
less than the 0.25 s bin time of our time-resolved spectra. 

Ejection of the accreted hydrogen likely requires that 
the accreted material first becomes highly stratified. In 
a well-mixed atmosphere, the H and He nuclei are effi- 
ciently coupled to the electrons (on which the radiation 
forces act) through Coulomb colisions, thus preventing 
any separation of the elements during the radius expan- 
sion episode. Variation in the H-fraction with depth may 
arise via a number of mechanisms, such as steady H burn- 
ing between the bursts or via convective mixing of deeper 
material during the previous X-ray bursts. However, as 
we show below, the properties of the X-ray bursts from 
4U 1636—536 place strong constraints on these mecha- 
nisms. 

In order to make a rough estimate for the requirements 
of the ejection of the H layer during a burst, we assume 
that all the radiation (emitted at ~ iEdd.n) during the 
radius-expansion episode is imparted as kinetic energy 
to the H layer. Since the H layer needs to be ejected 
in a short time te S 1 s, the total energy available is 
simply tgLEdd.H— 2 x 10'^^ ergs. This amount of energy 
can unbind a layer of pure hydrogen of total mass 

^ LEddMteRNS 

GA/ns 

_ -|^„18 / ^Edd,h^e \ / RnS \ 

\2x 1038 ergs / \10^ cm J 
( Mm Y' 

which extends down to a column depth of 

,j^6 f -^Edd,h^e \ / _RnS_\ 

~ \2x 1038 ergs / ^ 10^ cm J 

X fJ^y' gcm-^ (4) 

For the burst to be subsequently limited by the pure-He 
Eddington limit, the material below this column depth 
must be practically hydrogen free. This degree of strati- 
fication is rather difficult to achieve on an accreting neu- 
tron star, as we discuss b elow. 

■Gumm ing & Bildsten |200C1|D estimate the column 
depth at which hydrogen runs out to be 

.d ^ 7 X 10« ( - ) (^) (^) g cm- , 

\ommEddJ \ZcNoJ vo.7iy ^ 

where m is the mass accretion rate, Zcno is the mass 
fraction of CNO nuclei, and Xq is the mass fraction of hy- 
drogen in the accreting material. For this column depth 
to be comparable to the requirement derived in equa- 
tion the accretion rate or the mass fraction of CNO 
nuclei in 4U 1636—536 have to be very different than 
what is normally assumed. Depletion of hydrogen below 
a colu mn depth of c± 10^ gr cm~^ is possible, accord- 
ing to lFushiki fc Larnbl l)1987f) . if compressional heating 
is taken into account, the local accretion rate is compa- 
rable to the Eddington limit, and the core temperature 
of the neutron star is ~ 10^ K. This is again a rather 
implausible combination, since a high rate of accretion 



in a persistent source, such as 4U 1636—53 6, is inconsis- 
tent w ith a cool neutron-star core (see, e.g.. lBrown et alJ 
IT998I) . 

Convective mixing of deeper material during the pre- 
vious X-ray burst also appears to be incapable of reduc- 
ing the abundance of hydrogen at columns larger than 
~ 10^ gr cm ~^. Indeed, the conv ective zone in the 
simulations by iWooslev et al.l l)2004j) reached up to col- 
umn depths that were at least two orders of magnitude 
larger than required. The properties of the bursts in 
4U 1636—536 thus strongly suggest that there is an ad- 
ditional source of stratification which acts to separate the 
accreted elements, and allow the majority of the accreted 
H to be ejected during the brightest bursts. 

In Paper A, we also explored in detail the possibility 
that variations in -Fpoak could arise through systematic 
instrumental or analysis-related biases. For example, we 
considered the effect of deviations of the true emission 
spectrum from perfect blackbodies. Implicit in our es- 
timation of the bolometric fluxes in equation Q is a 
correction to the flux measured in the PCA bandpass; 
this correction adds around 7% to the peak 2.5-20 keV 
PCA flux of radius expansion bursts. Should the emitted 
spectrum deviate significantly from a blackbody, equa- 
tion (01 will not give the correct bolometric flux, poten- 
tially giving rise to spurious bolometric flux variations. 
As with 4U 1728—34, variations in the flux contribution 
from outside the PCA band of many times the typical 
bolometric correction would be required to account for 
the 30% variation in the peak flux of the majority of the 
radius expansion bursts from 4U 1636—536. We consider 
this unlikely. It is even less likely that such effects could 
account for the factor of « 2 variation we observe in the 
peak fluxes of all the PRE bursts from 4U 1636—536. 
Similarly, while variations in the persistent (non-burst) 
emission could in principle give rise to spurious scatter 
in the peak burst fluxes, the variation would have to be 
several times the total persistent flux level in the 2.5- 
25 keV band to account for the peak flux variation we 
observe. 

Although we found no evidence for quasi-periodic vari- 
ation in the peak fluxes of bursts from 4U 1636—536, or 
correlations between the peak PRE burst flux and flu- 
ence, we cannot completely rule out reprocessing as a 
contributing factor to the 7.6% peak flux variation of 
the bright PRE bursts. The lack of detectable quasiperi- 
odic variations in the peak burst fluxes could result from 
inadequate sampling provided by the burst times, or pos- 
sibly the intrinsic variability is aperiodic. As for the lack 
of a correlation between and i^pcak, the burst proflles 
from 4U 1636—536 were much more variable than those 
of 4U 1728—34, so that other (possibly related) varia- 
tions between these properties could perhaps mask an 
underlying correlation resulting from reprocessing. 

An alternative explanation for the low peak flux PRE 
bursts is that they arise from ignition of material con- 
fined to some fraction of the neutron star. Some au- 
thors have suggested that such confinement may occur 
as a result of the influence of the neutron star's magnetic 
fleld on th e fuel layer, or as a consequence of rapid ro- 
tation fe.g. lSztaino et al1ll985D . While recent modclhng 
indicates that r otation may instead enhance spreading 
ijSpitkovskv et al...2002.'l . variations in the properties of 
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bursts with M from some sources sugges t uneven distri - 
bution of fuel over the neutron star fe.g. lBildstenll20n(il) . 
This possibility ma y also help to explain the observation 
of doub le- ( Sztaino etaLlirOSSI and even triple-peaked 
l|van Paradiis et al.i '1986'^ bursts from 4U 1636-536. We 
also observed four double-peaked bursts in the RXTE 
sample. Confinement of the burst fuel could result in 
weaker radius-expansion bursts if the local Eddington 
limit was reached, but only over a fraction of the neu- 
tron star surface. In that case, expansion of the photo- 
sphere would still occur, at much lower total luminos- 
ity. If the expansion of the burning region of the photo- 
sphere went to sufficiently large radii that it appeared 
isotropic to a distant observer, the result might be a 
burst which appeared similar to a normal radius expan- 
sion burst, but at a significantly lower peak flux. This 
explanation appears inconsistent with the observed radii 
in the tail of the faint PRE bursts, which are in excess of 
that of some of the brighter PRE bursts. It would also re- 
quire that the factor of 1.7 between the peak fluxes of the 
two groups of PRE bursts was a coincidence. We note 
that several other explanations have be en proposed to 
explain t he bursts wit h multiple peaks l|Fiiiimoto et al.l 
[1988; Mch a fc Zvlstral ir9921. 

Previous estimates of the distance to 4U 1636—536 at 
between 6-7 kpc relied on the identification of the peak 
flux of the PRE bursts as iEdd,Hc ( Ebisuzaki 1987), in 
addition to a gravitational redshift measured from ab - 
sorption features in the burst spectra ijWaki et alJll984l) . 
While absorption features have now been detected in 
burst spectra from other sources using more modern, 
high-resolution spectroscopic instruments ijCottam et al.l 



|2002() . it is widely thought that the earlier detections us- 
ing proportional counters were more likely to arise from 
instrumental effects. Without the absorption line mea- 
surement, the analysis of [Ebisuzaki (1987) likely can- 
not substantially constrain the distance to 4U 1636—536. 
Despite the increase in measurement precision for peak 
burst fluxes with RXTE, the dominant uncertainty for 
distance estimates to most bursting sources remains the 
atmospheric composition, i.e. the value of X in equa- 
tion With the detection of PRE bursts that appar- 
ently reach cither iEdd,H or iEdd.Hc, the compositional 
question in 4U 1636—536 is resolved, and we can make 
distance estimates using either group of bursts. On cor- 
recting the observed fluxes for the gravitational redshift 
at the surface of a 1.4M0 NS with i? = 10 km, the 
fainter PRE bursts (which we assume reach LEdd.n) lead 
to a distance estimate of 5.95 ± 0.12 kpc. The distance 
estimate assuming that the brighter PRE bursts reach 
iEdd,Hc instead is fully consistent with this value at 
6.0 ±0.5 kpc, where the uncertainty arises from the stan- 
dard deviation of the peak fluxes for this larger group. 
From the faintest of this group of bursts we can derive 
a maximum distance of 6.6 kpc, or more conservatively 
(for a 2Mq NS) 7.1 kpc. 
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